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Abstract 

Fine powders often tend to agglomerate due to van der Waals forces between the 
particles. These forces can be reduced significantly by covering the particles with 
nanoscaled adsorbates, as shown by recent experiments. In the present work a quan- 
titative statistical analysis of the effect of powder flow regulating nanomaterials on 
the adhesive forces in powders is given. Covering two spherical powder particles ran- 
domly with nanoadsorbates we compute the decrease of the mutual van der Waals 
force. The dependence of the force on the relative surface coverage obeys a scaling 
form which is independent of the used materials. The predictions by our simulations 
are compared to the experimental results. 
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1 Introduction 



Dry granular powders are basic materials of great importance to the phar- 
maceutical industry. Bulk powders tend to agglomerate due to adhesive forces 
between the powder particles. This characteristic causes problems in the manu- 
facturing process of drugs, where accurate dosing is essential [1]. Furthermore, 
insufficient flowability can severely harm production devices as well as prevent 
inhalable drugs from reaching the lung's alveoli. 
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Fig. 1. (left) Model of a covered powder particle of radius R = 11 /jm and nanoad- 
sorbates of r = 50 nm (relative surface coverage of 2% on the lower sphere). Above, 
a second uncovered particle is approaching. The inset displays a zoomed view of 
the contact field, (right) Microscopic photo of a corn starch particle covered with 
aggregates of the nanomaterial Aerosil® 300, taken from [2]. 

Pharmaceutical powders typically have particle diameters on the micrometer 
scale and contact distances between the particles on the angstrom scale. If 
they are dry and not highly electrified, van der Waals forces are the strongest 
interparticle forces, exceeding gravitational and Coulomb forces significantly. 
As van der Waals forces are short-range it is possible to effectively decrease 
these adhesive forces by increasing the surface roughnesses of the powder par- 
ticles by covering them with nanoparticles (cf. Fig. 1). This was verified by 
tensile strength experiments, in which two powder layers were separated and 
the required force was measured [2,3]. 

H.C. Hamaker calculated the van der Waals forces between macroscopic spher- 
ical bodies by integration over all molecular dipole-dipole interactions [4]. Ac- 
cordingly 

f - A RlR2 m 



is the attractive interparticle force between two spheres of radii R\ and _R 2 , 
where all material constants are combined into the so called Hamaker constant 
A. H is the distance between the surfaces of the two bodies. In order to avoid 
a diverging force upon contact, one usually assumes that H cannot become 
smaller than a so called contact distance D = 4x 10" 10 m [5,6,3]. 

K. Meyer theoretically described the decrease of the tensile strength with in- 
creasing surface coverage assuming a stable three-point contact, i.e., effectively 
three nanoparticles act as spacers between two powder particles and combine 
them [2]. In bulk powders, two particles firstly joined by a one- or two-point 
contact are subject to torque, and their bonding remains unstable until a 
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three- (or more-)point contact is reached. 

Meyer considered an idealized system in which the three contact adsorbates 
are fixed on the vertices of an equilateral triangle [2,3]. This model allows one 
to calculate analytically the distance of the two bulk powder particles against 
the side length of this triangle. From this distance the interparticle force can 
immediately be derived using Hamaker's sphere-sphere model. 

The result of this calculation describes only two states sharply separated by 
a critical surface coverage p c . For p < p c the powder particles touch and 
attract each other strongly, while otherwise the nanomaterial establishes a 
gap between them and the short-range adhesive force vanishes nearly entirely. 
Obviously Meyer's model of equidistant spacing of the nanoparticles ignores 
any stochastic aspects which are relevant in reality. 

It is the aim of this paper to give a more realistic theoretical description of 
the effect of flow regulating nanomaterials on the interparticle forces acting 
between the larger powder particles. We simulate random coverages of the 
surfaces of the host particles with nanoparticles and calculate the van der 
Waals forces averaged over many realizations of randomness. We use spher- 
ical host particles and nanoadsorbates as a first approximation of reality, as 
suggested by microscopic photos of corn starch powder particles mixed with 
nanomaterials, see Fig. 1. Our generated data indicates a scale invariant be- 
havior, where an incomplete gamma function fits the resulting scaling function 
perfectly. Our results are compared to experimental data, demonstrating the 
improvements compared to Meyer's model. 



2 Simulation 



The simulation computes the averaged van der Waals force -F v dw versus the 
relative surface coverage p. The algorithm consists of three parts: (i) the prepa- 
ration of a host particle by covering it randomly with adsorbates, (ii) the po- 
sitioning of two host particles so that a stable three-point contact is formed, 
and (iii) the computation of the van der Waals force by summation over all 
resulting sphere-sphere interactions according to the Hamaker model. So each 
run of this algorithm gives one sample of the adhesive force which can occur. 

In the following the three steps of the simulation are shortly explained (for 
more detailed information see Appendix A). 
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Fig. 2. Algorithm for finding the three contact points, (a) Top view on the cov- 
ered host particle (white background) and (b) side view. The large dots in (a) are 
randomly distributed adsorbate particles. The uncovered powder (host) particle is 
approaching from above along the z-axis. The first adsorbate to be hit by the ap- 
proaching sphere is the one which is located nearest to the z-axis (Ci). Upon contact 
the circle in (a) describes the line of equal distance d in z-direction between the two 
host particles. C\ withstands the approaching sphere, so that the momentum of this 
moving sphere results in a torque (black arrow), (c) Due to the torque the upper 
host particle rolls off so that its center projected onto the x-y-plane moves along 
C±0, until this host sphere hits the second contact adsorbate (C2). (d) Now C\ as 
well as C2 hold the approaching host sphere so that its center (projected onto the 
x-y-plane) moves along the perpendicular bisector of the side C1C2, until a stable 
three-point contact is formed. 

2. 1 Preparation 

Initially, a host sphere is randomly and successively covered with nanopar- 
ticles, so that a uniform (Poisson) distribution with respect to the spherical 
surface area is established. Unlike the random sequential adsorption (RSA) 
model, our model does allow the adsorbate particles to overlap. Yet for the 
low coverages which we are interested in (p has the order of magnitude of few 
percent), the probability of two overlapping nanoparticles is negligible. 



2.2 Positioning 

In the second part of the algorithm a covered host and (for simplicity reasons) 
an uncovered powder sphere are brought to collide, cf. Fig. 1 left. The un- 
covered large powder sphere rolls off over the small adsorbates, until a stable 
three-point contact builds up (see Fig. 2); or, if the three potential contact 
candidates are too far away from each other, it rolls off until the two powder 
spheres directly touch each other. Knowing the characteristics of the contact, 
it is easy to calculate the new position s of the center of the uncovered powder 
sphere relative to the center of the covered host sphere, which is assumed to 
be located at the origin. 
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2.3 Computation of the force 



The coordinates Si describe the center positions of the Na adsorbates on the 
covered host particle. Let R be the radius of the powder particles, and r the 
common radius of all nanoadsorbates. A HH and Aah are the host-host and 
the adsorbate-host Hamaker constants, respectively. By vectorial summation 
over all sphere-sphere interactions according to Eq. (1) with two (host) powder 
particles, we get 



-FvdW — 
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(2) 



as the resulting interparticle force. 

Note that not only the three contact adsorbates but also all other adsorbate 
particles without direct contact to the uncovered powder sphere are taken into 
account. Most of them give only minor contributions to the adhesive force, 
as they are far away from the three-point contact and are therefore distant 
from the uncovered powder sphere. However, statistically some of them may 
come very close to the three-point contact and give force contributions of 
significant order of magnitude. Thus for a complete treatment of randomness 
it is reasonable to consider the influence of all the adsorbate particles for the 
computation of the force. 



3 Results and discussion 



For the further discussion, A = Ahh — Aah is assumed; most condensed 
phases, solid as well as liquid ones, have Hamaker constants of about 1CT 19 J 
interacting across vacuum, and similar values in air [7,6,3]. 

Concerning the adhesive effect between two particles, we have to examine only 
the absolute value of the van der Waals force, F = (|-F v dw|) as the mean of 
many samples generated by our algorithm. Some results of our simulations 
are shown in Fig. 3 (for details about the comparison with Meyer's model cf. 
Appendix B). 
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Fig. 3. Simulated data of the van der Waals force versus the relative surface coverage 
(A = 8 x 10~ 19 J, mean of 10 4 samples per point and standard error plotted). For 
comparison Meyer's model for R = 11pm and r = 50 nm is displayed. 

3. 1 Scaling 



3.1.1 Scaling properties 

For p — > 0, the force approaches the adhesion between two uncovered powder 
spheres, 

F - Fo = ^ oc R . (3) 



In Fig. 3 a linear right shift of the decrease of F in dependence of the surface 
coverage p for increasing r can be observed. This suggests the scaling form 



R \ r 



('-) ■ (4) 



The collapse of the data in Fig. 4 (top) shows that this scaling form works 
well for small values of pR/r, while for small F/R the curves become constant 
at different values, which grow roughly oc r/R. This seems plausible, because 
with a large gap between the two host particles the force is mostly affected by 
the three contact adsorbates only: 

min(F)^F m = 3-^^ 7 oc r (5) 



for r < R. 
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fi=5.5x10"V r=2.5x10" 8 m 
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R=2.2x10" 5 m, r=1.0x10~ 7 m 
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Fig. 4. (top) Semi-logarithmic plot of the rescaled simulated data (A = 8 x 1CP 19 J, 
mean of 10 4 samples per point). Increasing pR/r the curves saturate at a con- 
stant which scales approximately as oc r/R. The plot indicates a general behav- 
ior for pR/r — ► 0. (bottom) The converted quantity F = ^Zp of the same 
data produces a convincing data collapse. The inset and the right-hand plot indi- 
cate a power law behavior for very small pR/r and an exponential decay shortly 
before the minimum is reached, respectively. The scaling function is fitted by 
-0.499692 • (2 - T(3, 4.144 pR/r)). 

Our previous thoughts call for an improved ansatz of the form 



F = — 



F (R)-F(p,R,r) 
F (R)-F m (r) 



•(f) 



(6) 



where g is a scaling function in the range < pR/r < 2.5. Here the difference 
of the offset F to the data F is determined and divided by the possible 
maximum of this difference, given by F — F m . Scaling this way would give 
positive values with a maximum of 1. In order to obtain a function whose 
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curve represents the decrease of the adhesive force with increasing p in a more 
intuitive manner, the minus sign is inserted. Now, as can be seen in Fig. 4 
(bottom), we obtain a perfect collapse of the data. Note that F does not 
depend on A, i.e., it is independent of the used materials. 



3.1.2 Asymptotic properties of the scaling function 

For pR/r — > we expect a decrease of F which is oc p 3 , as can be explained 
as follows. In the limit of low densities, p — > 0, a decrease of the adhesive 
force can only occur if three potential contact adsorbates lie within a critical 
area between the two approaching host spheres. Assuming the positions of 
the adsorbate particles to be uncorrelated, the probability of this event grows 
oc p 3 . A logarithmic plot verifies this power law behavior, cf. the inset of Fig. 4 
(bottom). 

Increasing pR/r the semi- logarithmic plot in the right panel of Fig. 4 indicates 
a usual exponential decay until the minimum in F is reached at pR/r 2.5. 



3.1.3 Fit of the scaling function 

In our simulations the nanoscaled adsorbates are distributed randomly on 
the spherical host particle, i.e., they satisfy Poisson distributions in two di- 
mensions, projected onto a spherical geometry. In the analysis of 2D Poisson 
Voronoi cells it has been demonstrated that the distribution of cell areas can 
be fitted to gamma distribution functions with one, two or three parameters 
[8,9,10,11]. This suggests an ansatz of the form 

"(f)- (f)-(-f) 



for the derivative of the scaling function g, where a,/3 and 7 are constants. 
Integration yields 



pR 

pR 



g ^- — ^ = j a x 13 exp(— 72;) dx 




= «7^ 1 (r(/? + i)-r(/3 + i, 7 ^)) . (8) 

Here T(z) is the (complete) gamma function, and F(a, z) is the (upper) in- 
complete gamma function. 
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Fig. 5. Scheme of the tensile tester. 

Since F decreases as p 3 if p is small (see above), we expect the exponent (3 = 2. 
Indeed fitting Eq. (8) to our data in the range < pR/r < 2.5, one obtains 
(3 = 2.00 ± 0.01. Setting (3 = 2 and performing another fit, we determined the 
fit parameters as a = —35.56 ± 0.07 and 7 = 4.144 ± 0.003. So the numerical 
result for the scaling function is 



=-0.500. ( 2 -P (3. 4.14 £)) 



(9) 



In Fig. 4 (bottom) one can see that it matches the rescaled numerical data 
perfectly. 



3.2 Comparison to experiments 



The theoretical results can be compared to previous experiments with a tensile 
tester [2], cf. Fig. 5. 

In this experiments a vaselinated measuring plate is brought into contact with 
the surface of a powder sample. Lifting the plate a tensile force is measured 
which allows one to determine the interparticle forces needed to separate two 
powder layers (for details see Appendix C). 

A comparison of the experimental data with the results of the simulation is 
shown in Fig. 6. Here the quantity of the simulated adhesive force F is arbitrar- 
ily scaled against the quantity of the experimentally measured tensile strength 
a so that the offsets collapse. For p = there are no nanoscaled adsorbates 
on the host particle, and so F and a are generated only by powder-powder 
interactions. Although the experimentally measured forces for p > are sys- 
tematically larger than the predicted ones, they show a qualitatively similar 
type of curvature, confirming the approach of studying random coverages. 

At present the origin of the systematic deviations is not yet entirely clear. 
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Fig. 6. Qualitative comparison of the simulated adhesive force F with experiments in 
which the tensile strength a was measured; note that these are two different physical 
quantities F and a, scaled against each other arbitrarily in order to give combina- 
tions showing common trends in simulation and experiment. The experimental data 
is taken from [12]. (top) After some mixing, by which agglomerates of the nanoma- 
terial are crushed into smaller aggregates, aggregates of Aerosil® 300 deposited on 
corn starch powder particles (R = 11 /xm) have average radii of r = 42 nm. The sim- 
ulation used A = 6.5 x 10~ 20 J as Hamaker constant for this material combination, 
(bottom) Analogous graph for Aerosil® OX 50 on corn starch with R = 11 /im, 
r = 80 nm and A = Q.5 x 10~ 20 J. 

They could be attributed to the influence of additional liquid bridges between 
the powder particles. However, with corn starch the formation of liquid bridges 
is only observed at relative humidities of more than 40% [12], and all experi- 
mental measurements were performed well below this critical limit. 

Aggregation is another important aspect. The adsorbate particles mentioned 
above are aggregates of nanoscaled primary particles [2]. Depending on the 
number of combined primary particles the diameters of the adsorbates vary. 
The radii r used in our simulations are averages over many adsorbates whose 
sizes were measured by atomic force microscopy during the experiments. Ex- 
perimentally the measured radii have errors of the order of about r/2 [12]. 
However, in the simulation statistical fluctuations of the adsorbate radii are 
not included in order to keep the problem simple and fast computable. In 
this way we make systematic errors, but a relative error in r of about 0.5 
alone cannot explain the quantity of the deviations between simulation and 
experiment, cf. Fig. 6. 

The precise three-dimensional shapes of the adsorbates are also not yet com- 
pletely analyzed. Experimentally the aggregates are placed onto the surfaces 



10 



of the host particles by simply mixing the nanomaterial with the powder. The 
number of adsorbed aggregates increases during mixing, and the surface cov- 
erage of the host particles grows. As an unintended side effect of this mixing 
a 'grinding' of the aggregates occurs [1]. While at the beginning of the mixing 
process we expect the shape of the adsorbed aggregates to be spheric as an 
approximation, their shape gets continuously flatter due to the mechanical 
stress during mixing. So we expect a deviation upwards from the computed 
adhesive force because flat adsorbates generate a smaller distance between two 
powder particles than the spherical adsorbates assumed in the simulation. 

Additionally, it seems to be likely that the experimental method overestimates 
the actual forces. In fact, the top layer separated by the tensile tester is not a 
monolayer of spheres. Instead we expect a rough and irregular surface, espe- 
cially at the edges of the sample, and so the force needed to produce such a 
surface is presumably higher. 

Nevertheless the systematic nature of the deviations allows us to conclude that 
the simulation results may serve as lower bounds to the actual forces and to 
the minimum surface coverage needed to optimize the powder flow. 



4 Conclusions 



In the present work we have investigated how van der Waals forces between 
two spherical particles are reduced by randomly deposited nanoscaled adsor- 
bates. In contrast to an earlier work by Meyer, who considered a triangular 
coverage, the statistical analysis of randomly deposited nanoadsorbates leads 
to an improved prediction of the force reduction as a function of the coverage. 
More specifically, a continuous curve instead of a step-like function is obtained. 
Varying parameters we have identified scaling laws which are independent of 
the used materials. A comparison with previously recorded experimental data 
shows that the force predicted by the simulation may be used as a lower bound 
for the actual force, though the experimental data is not sufficient and not pre- 
cise enough to verify the functional behavior properly. This could be a starting 
point for further experiments. Future simulations could include a statistical 
distribution of the adsorbate size. However, the decline of the adhesive force 
computed in our model may act as a lower approximation of the minimum 
surface coverage which is needed to optimize the powder flow. 
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A Simulation in detail 



A . 1 Preparation 

For small p there is no difference - respecting the problem - whether to regard 
two covered powder particles, or one uncovered particle and one which is 
covered at double density, cf. Fig. 1. Thus, for simplicity, the preparation part 
of the simulation covers only one host particle stochastically with the coverage 
density 2p. 

The center of the sphere to be covered is the origin in our coordinate system. 
Random sphere point picking must be done with respect to the sphere's sur- 
face area. In terms of spherical coordinates, where 9 is the polar and ip the 
azimuthal angle, 



is an area element on the unit sphere. So a uniform point distribution is 
achieved by rectangular distributions of 



As a further simplification, we cover only one half of the first sphere, where 
the contact to the second, uncovered sphere will be made. The neglected ad- 
sorbates would hardly contribute to the bonding force, as they would be too 
far away from the second powder particle. 

In this way adsorbates are successively placed upon the host, with their centers 
at coordinates &i,(pi, until the desired coverage is reached. Let Os = 4-7T-R 2 be 
the host's surface area and Ca = ?rr 2 the area of perpendicular projection of 
an adsorbate on this surface. So the number of adsorbates on the host sphere 
can be calculated from the relative surface coverage p as 



dVt = dip d6 sin 9 = — dip d(cos 6) 



(A.l) 



pE[0,2n) and cos0 G [-1, 1] 



(A.2) 




(A.3) 



Note that the factor 2 due to covering only one host particle and the factor | 
due to covering only half this particle cancel each other out. 
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Fig. A.l. (left) If the circumcenter c of the contact adsorbates is more than D/2 
away from the surface of the covered powder particle (located at the origin O), the 
second powder particle is placed in such a way that c bisects the distance between 
these two powder particles, (right) If c is less than D/2 away from the surface 
of the covered powder particle, i.e., that the two powder particles will touch each 
other, the second powder particle is placed in the direction of c so that there is just 
the contact distance D between the two powder particles. 

A. 2 Positioning 

After complete preparation of the first sphere, the second uncovered sphere 
approaches along the z-axis in negative direction, until it hits the covered 
hemisphere. In order to determine the three contact points as easy as possible, 
the coordinates of the Na adsorbates are transformed by projection on the x- 
y-plane, (x,y,z) — > (x,y), which is a sufficient approximation for r <^ R. 

The roll-off mechanism described in Fig. 2 and its caption yields the three 
contact points; consequently the contact position of the second powder par- 
ticle can be calculated, which will be explained in the following. With the 
coordinates of the centers of the three contact adsorbates the circumcenter c 
of the triangle spanned by these three points can easily be constructed. Note 
that the center of the covered powder particle was defined to be located at 
the origin. So for symmetry reasons, see Fig. A.l, the center position s of the 
second, uncovered particle can be computed as 

2 c if \c\>R + § 

s=l ' . (A.4) 

(2R + D)±c if \^<R+§ 

As the algorithm for finding the three-point contact does not check for a direct 
contact between the two host particles, the artificial definition of s for |c| < 
R+ y must be introduced. So we ensure a minimum distance of D between the 
two host particles. For this case of contacted powder particles the adsorbates 
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Fig. B.l. Surface layout according to Meyer's model. The nanoadsorbates (large 
dots) occupy the vertices of the equilaterally triangulate surface of the powder 
particle (white background). Distributing the surface area equally to the adsorbates, 
the area of the highlighted cell belongs to each individual one. 

are of minor importance, since the direct powder-powder interaction exceeds 
any powder-adsorbate interaction by some orders of magnitude. 



B Conversion of the parameters from Meyer's model 

In order to be able to compare Meyer's model with our simulated data, one 
has to convert the control parameter used in this study, the relative surface 
coverage p, into the quantity x\ used by Meyer, cf. Fig. B.l and [3,2]. 

Eq. (A. 3) yields the total number of adsorbates on one host particle, p \— > Na- 
Although it is not possible to cover the surface of a sphere globally by a 
hexagonal lattice it is assumed that such a structure can be established locally. 
According to Meyer's model, adsorbates are set on the vertices of this lattice, 
so that a uniform pattern of equilateral triangles evolves. To each adsorbate 
belongs a cell of the area 



cf. Fig. B.l. With elementary geometrical considerations one can calculate the 
side length of an equilateral triangle of the area as 




(B.2) 



Consequently 



a — 2r 



a 



(B.3) 



x 1 = 



2 
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Nanoma/tericd 


Diameter of the primary particles 
[nm] 


ooiiQ density 
[g/cm 3 ] 


opecmc sunace 
[m 2 /g] 


Aerosil® 300 


7 


2.2 


300 ± 30 


Aerosil® OX 50 


40 


2.2 


50 ± 15 



Table C.l 

Specific properties of the used nanomaterials. 

follows, which can be put in the formulas for h(xi) and F v dw{h), see Eqs. (27) 
and (24) in [3]. Hereby a mapping p \— > F v dw gets feasible. 



C Experiments 

C.l Materials and methods 

Corn starch, Cerestar® GL 03406, was purchased from Cerestar Germany, 
Krefeld. The spherical particles have diameters of 22 pm. They are very co- 
hesive. As nanomaterials we used Aerosil® 300 and Aerosil® OX 50 both 
supplied by DEGUSSA AG, Hanau, Germany. Their specific properties are 
summarized in Table C.l. 

Both types of Aerosil® are hydrophilic. Due to their small sizes they highly 
agglomerate. 

In our experiments we determined the tensile strengths of powder mixtures 
consisting of corn starch and 0.2% of a nanomaterial. In order to prepare these 
mixtures first the nanomaterial was given into a glass vial with a volume of 
500 ml. Then in a second step 100 g of corn starch were added. The degree of 
filling was given by 0.4. The powders were mixed by means of a Turbula mixer 
(Type T2C, Nr. 950353, W. Bachofen AG, Basel, Switzerland) at 42rpm. A 
comparison of the particle sizes of the mixture components suggests that the 
mixing process involves a grinding of the agglomerates formed by the nano- 
materials. 

C. 2 Measurement of tensile strengths 

The interparticle forces acting in the mixtures consisting of corn starch and 
0.2% of a nanomaterial were measured by means of tensile strength tester 
developed by Schweiger [13,14] and modified by Anstett [15], Fig. 5. 

The powder to be tested is filled into the sample container having a volume 
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of 4.62 cm 3 (diameter: 3.43 cm, height: 0.50 cm). By means of a sieve with a 
sieve size of 315 /im it is sieved directly into the sample container. In order 
to achieve a homogenous filling the filled container is gently tapped onto a 
soft pad. By scraping off the surplus testing material by means of a slide a 
perfectly flat surface of the powder bed can be obtained. This step has to be 
performed very carefully in order to avoid even the slightest compression of 
the powder. 

The sample container prepared in this way is placed on a small table which can 
be moved up and down by means of the linear positioner M-410-21 (Physik 
Instrumente GmbH&Co, Waldhorn). By coarse adjustment screws the whole 
unit consisting of a table as well as a linear positioner is moved upwards so 
that only a very small gap of about 1 to 2 mm remains between the powder 
surface and the measuring plate. From this position on the further movement 
of the table by the linear positioner is controlled by the controller C-832.00 
(Physik Instrumente GmbH&Co, Waldhorn) which runs under the software 
package Pro Move v2.20. The steps performed by the linear positioner can be 
simultaneously seen on the screen of the PC. With a vertical speed of 3.4 fim/s 
the table is moved upwards until the measuring plate comes into contact with 
the powder surface. The movement is stopped as soon as a defined negative 
weight of the measuring plate is measured. 

The cylindrical measuring plate has a lower surface of 0.385 cm 2 . By means 
of a frame made from Aluminum it is fixed on a load cell ZER 10 (Wipotec 
GmbH, Kaiserlautern, Germany). The measuring range of this load cell covers 
a range of 2g with a resolution of 0.1 mg. The force acting on the load cell 
is registered by the software Wipotec- Terminal v3.78 at a scan rate of three 
measurements per second. The signals of the load cell are mass calibrated. In 
order to obtain the corresponding forces the signals have to be multiplied by 
the gravitational constant g. 

Before each measurement a solution of 5% (w/w) white petrolatum in petro- 
lether is sprayed onto the lower side of the measuring plate. After evaporation 
of the solvent the measuring plate is coated with a thin layer of petrolatum 
as a sticking agent to ensure that the powder particles adhere on its surface. 

After a pause of 10 s the table with the sample container is moved downwards 
at a speed of 1.7 m/s. Due to the interaction between the layer of powder 
particles fixed at the lower side of the measuring plate and the particles re- 
maining in the adjacent powder layer in the sample container the load cell 
measures a force. This force increases with continuing downwards movement 
of the sample. It reaches its maximum when the two powder layers separate. 
The force measured by the load cell drops to that force corresponding to the 
weight of the powder layer adhering at the lower surface of the measuring 
plate, Fig. C.l. 
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Fig. C.l. Time course of the force during a measuring cycle. 



In order to calculate the tensile strength a the force Fp measured at point F 
is subtracted from the maximum force Fp measured at point E. This value is 
then divided by the lower surface Am of the measuring plate: 
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